Numerical study of the Coulomb blockade in an open quantum dot 
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The Coulomb blockade in an open quantum dot connected to a bulk lead by a single mode point 
contact is studied numerically using the path-integral Monte Carlo method. The Coulomb oscillation 
of the average charge and capacitance of the dot is investigated, and is compared with the analytic 
expression for strong tunneling. At the degeneracy point, we observe logarithmic divergence of the 
capacitance for strong backscattering at the point contact. This observation supports the conjecture 
that the nature of the present system at the degeneracy point is described by the two-channel Kondo 
problem for an arbitrary strength of tunneling. 
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In a mesoscopic structure called a quantum dot, in- 
teraction between electrons restrict the amount of the 
inner charge Q at low temperatures T ^ Ec, where 
Ec = e^/2Co is the charging energy, and Co is the static 
capacitance of the dot. The electrostatic energy of the 
dot is given as Eq — (Q — eN)^/2CQ, where N is propor- 
tional to the gate voltage Vg. If transport of electrons is 
governed by weak electron tunneling between leads and 
a dot, one observes sharp periodic peaks of conductance 
as a function of the gate voltage. Here, the peaks corre- 
spond to the degeneracy points = n+1/2 between two 
charge states, Q = ne and Q = {n+ l)e. Away from the 
peaks, the amount of the inner charge is fixed as Q = ne, 
where n is an integer minimizing Eq. These phenomena 
are called the Coulomb blockade. For strong tunneling, 
which is realized in a quantum dot made in semiconduc- 
tor heterostructures, the Coulomb blockage is suppressed 
by large charge fluctuation in the dot. 

In order to deal with the strong tunneling regime, it 
is convenient to consider a system consisting of a dot 
and a half-infinite lead connected by a single mode point 
contact. The Coulomb blockade in such a structure can 
be observed by measuring the average charge and the ca- 
pacitance, and have been studied both theoreticallyii^i^ 
and experimentallyi^ The peculiar nature of this system 
is logarithmic divergence of the capacitance character- 
ized by the analogy to the two-channel Kondo problem®; 
Matveev has revealed that in the present system loga- 
rithmic divergence of the capacitance appears at the de- 
generacy points in the regime of both weak and strong 
tunneling through the point contact The consistency 
in these opposite limits suggests that the Coulomb block- 
ade system at the degeneracy points can be described 
effectively by the two-channel Kondo model for an arbi- 
trary strength of tunneling. One way of confirming this 
conjecture is to numerically study the intermediate tun- 
neling regime. Recently, the numerical renormalization 
group (NRG) method has been applied to this problem,^ 
and nonmonotonic growth of the Coulomb blockade has 
been reported in the strong tunneling regime. However, 
this calculation is based on a model of the weak tunnel- 
ing regime. It is quite necessary to study this problem 
by using a model applicable to both weak and strong 



tunneling. 

In this paper, we study the Coulomb blockade phe- 
nomena for an intermediate tunneling in an open quan- 
tum dot — a dot strongly connected to a bulk lead. If 
the curvature of the constriction near the center of the 
point contact is smooth, the system is essentially one- 
dimensional, and described by the Tomonaga-Luttinger 
liquid (TLL): an effective model of low energy excita- 
tions in an interacting one-dimensional electron system. 
According to Ref. [3, the Hamiltonian of the system is 
bosonized as 
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where we assume that the dot (lead) corresponds to the 
region a; > (a; < 0). The subscripts p and a stand for 
a charge and spin mode, respectively. The first term in 
r.h.s. of Eq. ([T]) describes a bulk spinful TLL with the 
sound velocity . The positive parameter represents 
an interaction strength: Ki, < 1 for the repulsive case; 
Ki, > 1 for the attractive case; K^, = 1 for the noninter- 
acting case. Charging effect in the dot and backscattering 
of electrons at the point contact correspond to the other 
two terms in Eq. ([1]), whose coefficients are defined as 
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Here r is the reflection amplitude, D = 2tt/At is the 
bandwidth cutoff, and N = CoVg/e is the dimensionless 
gate voltage. We note that the present model can cover 
the whole region between the weak tunneling limit (r 
0) and the strong tunneling limit (r — )■ 1). 

Since the charging and backscattering term in Eq. ([1]) 
are expressed only by the fields at the origin, we can 
integrate out the boson fields in the bulk part and obtain 
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the effective action 

S = So + Sc + Sv, 



EE 



Sc = UJdT [0p(t) - ^iV]2, 
Sy = V / dr cos (/)p (r) cos (t) , 



(3) 
(4) 

(5) 
(6) 



where the spatial coordinate x is omitted, and 4>„ denotes 
the Fourier component of </)^. The path- integral Monte 
Carlo (PIMC) method^"^- is one of the powerful ways to 
simulate quantum systems described by an effective ac- 
tion such as Eq. ^ . Recent development^ of this method 
has remarkably improved the efficiency of the simulation 
at low temperatures. Our Monte Carlo simulation is con- 
structed by local update in the Fourier space and global 
update based on the cluster algorithmi^i^ Similar simu- 
lation has been performed in our previous work,^*' which 
has discussed the impurity problem in a spinful TLL cor- 
responding to the case of C/ = 0. 

Discretizing the imaginary time into L steps, we define 
the jth step on a path as 4>,yj = 0^(j/3/L), and its Fourier 



transform as 
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thus obtain the discretized effective action suitable to the 
local update as 
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where Ar = In a local update, for each pair of 

v and k a value of (pi^k is randomly chosen following a 
normal distribution cx e^l' 
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(fc^0,L/2). (11) 



Then a new path is obtained from the inverse 
Fourier transform, and is accepted with probability p = 
minllje""^"^*^}. The local update in the Fourier space 
globally changes the paths in the real space. In the pres- 
ence of strong backscattering, however, relevant update 
tends to be rejected at lower temperatures, because such 
global changes do not consider the potential configura- 
tion. 

The basic idea of the cluster algorithm is as follows. A 
reflection mirror is suitably located, and every two sites in 
a discrete path are connected with a certain probability. 



Then reflection of the connected sites against the mirror 
is accepted with a probability determined by the cost of 
reflection of the cluster. If we rewrite the effective action 



as 



Sc = UATY}<l>Pj-^Nf, 
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with the kernel defined by 



k=-L/2+l 



L/2 

iyt|g(2^Vi)(i-i')fe^ 



(15) 



we can design two types of cluster update as follows. In 
a smgZe-field cluster update, clusters for the charge and 
spin field are constructed separately with the bond prob- 
ability between the jth and j'th site 



Puj]' = max{0, 1 - c-xp{~2K^jjiip^jip^f)}. 



(16) 



Here we define the relative field Lp^j = (p^j — ■nM^ mea- 
sured from the mirror located at (j)^ = ttMi, with an in- 
teger Mi^. In a double-Held cluster update, on the other 
hand, a pair of the jth and j'th site for both the charge 
field and the spin field is connected using the probability 



Pjj> =max<^ 0, 1 - expl -2^^Ki,jj>Lp^jipi,j 



(17) 



with a half integer M^. Clearly, Sy remains unchanged 
after these global moves, and the strength of backscat- 
tering does not matter. Sc is invariant exceptionally 
in the single-field case for the spin mode, but it gener- 
ally changes in the other cases. Then the cluster up- 
dates for the latter cases are accepted with the probabil- 
ity p = min{l, e"'^'^'^}, where the change in Sc has the 
form 



ASc = AnUAriN - Mp) 
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The kinks, which run from a potential minimum to an 
adjacent one of the double-cosine potential in Eq. (|14p . 
are inserted in the path by these two types of cluster 
updates. These kinks are relevant structures when one 
considers strong backscattering (corresponding to weak 
tunneling) at low temperatures. 

In what follows, the results of our PIMC simulation 
are presented. We use [/ = 1 as the energy unit, and fix 
the width of the time step to Ar = 1/4. A noninteract- 
ing TLL {Kp = Kr, = 1) is assumed in the bulk part. 
Because of the periodicity of the double-cosine potential 
in Sv [see Eq. ([T^ ]. we can restrict the value of N in 
the range < < 1. Since the charge in the dot is 
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bosonized as Q = {G/iT)(f>p{i 
measured as 

(Q) = -(0p), 4>p 



0), the average charge is 



1 



dr (t>p{T), 



(19) 



where (O) = Z^^ J V/ppVcji^Oe^^ denotes the expecta- 
tion value of an observable O. In Fig. [U we plot (Q) 
for y = 4 and 8 as a function of the gate voltage N. 
As the temperature decreases (the inverse temperature 
L increases), the Coulomb blockade is enhanced, and the 
maximum and minimum point of the values of (Q) — eN 
approaches = 1/2 for each V. If we focus on the low- 
est temperature L — 400, the curve appears to converge 
except near the degeneracy point N = 1/2. Therefore, 
we can expect that the data &i L — 400 is comparable to 
the analytic result for T = in the strong tunneling limit 
{y ^ Q). In Ref. 2, the charge field (/)p is pinned down 
to ttN due to the large charging energy T ^ Ec < D, 
and the backscattering strength V is renormalized by the 
charge fluctuation. As a result, the average charge is cal- 
culated for small V (corresponding to strong tunneling) 
as 



(Q) 



= N - 



log(e |r p cos^ irN) sin 27rA^, (20) 



where 7 = e^ with the Euler's constant C ~ 0.5772. The 
inset of the Fig. [T] shows the comparison for V — 4 and 8 
between our data of L = 400 and Eq. ^U\j. For V ^ A, the 
analytic result shown as a solid line is almost superposed 
on our data points; the expression ((20|) based on the pin- 
ning of (j)p describes well the Coulomb oscillation of (Q). 
For V ^ 8, on the other hand, our data points away from 
N — 1/2 disagree with the analytic expression; the sys- 
tem is no longer in the strong tunneling regime for such 
a large value of V. It should be noted that, the Coulomb 
blockade grows monotonically with increasing V and L, 
and that the reentrant behaviors in the strong tunneling 
regime reported in the NRG results^ are not observed in 
our simulation. 

Near A^ = 1/2 the oscillation curve in Fig. [T] slightly 
depends on L at lower temperatures. This is related to 
the singularity derived from the two-channel Kondo prob- 
lem hidden in the Coulomb blockade system at the de- 
generacy point. To see the peculiar nature of the system 
more clearly, we now measure the capacitance of the dot 
defined as 
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(21) 



In Fig. [2I we plot the capacitance ior V — 4 and 8 as 
a function of the gate voltage A^ at low temperatures. 
The amplitude of the Coulomb oscillation of the capaci- 
tance grows larger as L and V increase. The temperature 
dependence strongly depends on the gate voltage; away 
from A^ = 1/2, the capacitance is convergent at low tem- 
peratures whereas the peak at A^ = 1/2 keeps growing 
and sharpening up to the lowest temperature L = 400. 
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FIG. 1: (Color online) Coulomb oscillation of the average 
charge of the quantum dot for different temperatures L — 10, 
20, 50, 100, 200, and 400. A bundle of curves with smaller 
amplitude corresponds to 1^ = 4, and the other to 1^ = 8. 
Inset; Comparison of the average charge for V = 4 and 8 
between our PIMC data at L = 400 (crosses) and the analytic 
result for small V at zero temperature (solid line). 



3 
2.5 

2 
1.5 

1 

0.5 





FIG. 2: (Color online) Coulomb oscillation of the capacitance 
of the dot for different temperatures L = 10, 20, 50, 100, 
200, and 400. A bundle of curves with smaller amplitude 
corresponds to V = 4, and the other to V = 8. 



In Fig. [21 we show the capacitance for intermediate 
tunneling = 8 as a function of the inverse temperature 
L for different values of the gate voltage A^. In this figure, 
the temperature dependence of the capacitance is clas- 
sified into three different types: when the temperature 
decreases, (I) for A < 0.40, the capacitance decreases 
monotonically; (II) for A^ = 0.45, the capacitance once 
increases, and decreases below a crossover temperature; 
(III) for A^ = 0.50, the capacitance increases monotoni- 
cally. In Ref. @, it is suggested from analytic results that 
the low-energy property of the Coulomb blockade system 
near the degeneracy point (A^ = 1/2) is governed by the 
physics of the two-channel Kondo problem for an arbi- 
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FIG. 3: (Color online) Temperature dependence of the capac- 
itance of the dot for V — 8. The gate voltage N is varied from 
0.15 to 0.50 (from bottom to top) at intervals of 0.05. 
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FIG. 4: (Color online) Logarithmic divergence of the capaci- 
tance of the dot for A*' = 1/2. The capacitance is plotted as a 
function of the inverse temperature L. The logarithmic fit of 
the data points for each value of V is also shown by a dotted 
line. Inset: l/-depenence of the coefficient of the logarithmic 
divergence. The solid line shows the analytic result. 



trary strength of tunneling. This prediction is supported 
by the logarithmic divergence of the capacitance for the 
type III. The transient increase in the capacitance for the 
type II is also a sign of the two-channel Kondo physics. 
For N ^ 1/2, an intrinsic low energy cutoff in the strong 
tunneling limit (V 0) is given bj^ 



(22) 



characteristic of the two-channel Kondo model disap- 
pears. One can confirm in Fig. [3] that the crossover 
temperature in the type II is roughly given by F. 

In order to study in detail the logarithmic divergence 
of the capacitance, we focus on the degeneracy point 
N — 1/2 and show the capacitance in Fig. [4] as a func- 
tion of the inverse temperature L. The results for dif- 
ferent backscattering strengths up to F = 8 are shown. 
One can see that our data points diverge logarithmically 
with decreasing temperature for any value of V as far as 
our simulation is performed. Finally, we compare this 
logarithmic behavior of our data with the analytic result 
obtained in the strong tunneling limit (T^ — > 0). Accord- 
ing to the analytic expression,^ the capacitance diverges 
at low temperatures as 



4^ 

TT ^\ T 



(23) 



At low temperatures below F, the logarithmic divergence 



We can directly compare the coefficient of the logarithm 
between Fig. Hand Eq. In the inset of Fig. IH the 

coefficient derived from the logarithmic fit is plotted as a 
function of V. For small V, i.e., in the strong tunneling 
regime, our data points agree well with the analytic result 
47r7|r|^Co/e^ ~ 0. 0044291^^ shown as a solid line. For 
larger values of V, however, they grows faster than the 
analytic result. Thus, the analytic expression ([23| fails 
for intermediate strength of tunneling. Nevertheless, the 
logarithmic divergence in Fig. 2] in this regime clearly 
indicates that the present system is understood by the 
physics of the two-channel Kondo model for any value of 
V. 

In summary, we have started with the bosonized 
Hamiltonian of an open quantum dot, and demonstrated 
the PIMC simulation of the Coulomb blockade phenom- 
ena at low temperatures. The Coulomb oscillation of the 
average charge and capacitance of the dot is studied in 
the intermediate tunneling regime, where the analytical 
result for strong tunneling is invalid. The amplitude of 
oscillation monotonically increases with decreasing tem- 
perature or increasing the backscattering strength. We 
confirmed that logarithmic divergence of the capacitance 
appears at the degeneracy point N = 1/2 even in the in- 
termediate tunneling regime. This observation indicates 
that the present system has a two-channel Kondo nature 
for an arbitary strength of tunneling. 
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The variance a^o = oo means that a uniform distribution, 
e.g., ranged from — tt to tt is used for a local update of (pao- 



